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We investigate interaction effects in transport phenomena in bilayer graphene (BLG). For the 

minimal conductivity in pristine BLG, wc find that the conductivity assumes a constant value in 
the limit T — > 0, with the first correction being oc \/T. This has to be contrasted from the standard 
1/T^ in Fc :nni liquids (neglecting additional logarithms and above all disorder). We furthermore 
study the Coulomb drag resistivity between two BLGs in the whole range from deep within the 
Fermi liquid regime all the way to the charge neutrality (CN) point. We find that in the Fermi 
liquid regime drag behaves very similarly to drag in a standard two-dimensional electron gas. In 
contrast to monolayer graphene, we find no saturation of drag as a function of the distance d for 
realistic parameters. In the vicinity of CN, wc find an interesting interplay between interaction 
effects and disorder, like in the case of monolayer graphene. Here the drag resistivity strongly 
depends upon the ratio of the corresponding scattering times. 



I. INTRODUCTION 

It is well know that decreasing dimensionality is a 
means of increasing the effect of interactions in electronic 
systems. In one dimension this leads to Luttinger liquid 
behaviour and a variety of manybody instabilities.^ In 
two dimensions there also is a huge variety of strongly 
interacting electronic systems,^ leading to phenomena as 
diverse as the fractional quantum Hall effect^'^ or high 
temperature superconductivity^'®, to name two of the 
more prominent ones. A purely two-dimensional system 
by construction is monolayer graphene (MLG),^'* which 
is a hexagonal arrangement of carbon atoms. This sys- 
tem has remarkable electronic properties,^ and a lot of 
these can be traced back to the fact that the low-energy 
theory assumes the form of an emergent massless Dirac 
equation. In this system, while with magnetic field the 
fractional quantum Hall effect has been observed, -'^'^ with- 
out magnetic field the observation of interaction effects 
is scarce to date. One of the reasons of the robustness 
of the underlying electronic system against interactions 
is the vanishing density of states at the Fermi level. 

In bilayer graphene (BLG) systems the situation is 
different. ^'^ The major difference stems from the fact that 
instead of massless Dirac particles, we deal with massive 
quasi-particles."'"* These occur since in the low-energy ef- 
fective theory we have two bands which touch at one 
point, however not in a linear but in a quadratic fash- 
ion, implying a finite band mass."'^^ As a consequence, 
even at charge neutrality (CN), the density of states is 
finite instead of zero as in MLG. This implies that the 
system is much more unstable with respect to interac- 
tions. A variety of different symmetry broken states with 
a finite gap are conceivable,^^, for instance by explicitly 
breaking the interlayer symmetry by an external gate 
voltage. ^®'^'' Without external perturbations there are in- 
dications that below T = 5 K BLG undergoes a sponta- 
neous transition to a gapped state. The exact nature 
of this state is currently under dispute^^"^^ and further 
theoretical and experimental activities are required to 
identify it unambiguously. 



Within this paper we do not attempt to speculate 
about the nature of this symmetry broken state but in- 
stead focus on interaction effects and their signature in 
transport properties. Consequently our theory only ap- 
plies for temperatures which are above the gap scale. It 
also neglects the scattering of electronic degrees from col- 
lective modes which can be precursors to the transition 
to the ordered state. 

In a first part, we investigate the minimal conduc- 
tivity of clean BLG. As in the case of MLG there is 
a well-defined minimal conductivity, even in absence 
of impurities solely due to (screened) Coulomb interac- 
tion. This is a consequence of the particle-hole symme- 
try in the model and holds true for every particle-hole 
symmetric system. We study the minimal conductivity 
as a function of temperature, and find that it follows 

Cmin = X (^'''■^ + ^'^■^^V^) ^^"^ reasonable tem- 
peratures T. In this expression m is the effective mass 

of the BLG quasi-paticles, and e is the dielectric con- 
stant of the environment. Interestingly, it saturates for 
low temperatures as opposed to the standard behavior 
in Fermi liquids where it diverges (due to the vanishing 
phase space). 

In a second part we study the effect of Coulomb drag. 

This is a very direct measurement of Coulomb interaction 
between two electrically isolated two-dimensional elec- 
tronic systems. The experimental setup goes as follows: 
A current is driven in one of the layers, called the active 
one. Then Coulomb interaction can transfer momentum 
to the; other layer, called the passive one. If the cur- 
rent is not admitted to flow there, a voltage drop will 
be induced. This can be measured, and the ratio of the 
current in the active layer and the voltage drop in the 
passive layer is called the drag resistance. 

Our study provides the first comprehensive study of 
Coidomb drag in BLG in all limits ranging from CN to 
the Fermi liquid (FL) regime. To the best of our knowl- 
edge experimental results for such a setup are currently 

unavailable, and we are only aware of one theoretical 
study.^"^ For high doping we find that the system be- 
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haves like a standard FL, and we discuss all the limiting 
cases as a function of temperature T and distance d be- 
tween the layers. Unlike in the case of MLG, where there 
is a well-defined regime of 'zero' distance, this regime is 
lacking here due to the small screening length. In the 
opposite limit, close to CN, we find that drag depends 
very sensitively on the ratio between relaxation due to 
interaction and disorder, thus inelastic and elastic scat- 
tering. A similar behavior has recently been observed in 
theoretical studies of drag in parallel MLGs.^^'^^ 

Technically, this paper follows in large part the lines 
of a recent publication^'' on drag in a MLG setup and 
consequently we present little technical detail. The gen- 
eralization to the present situation is pointed out in the 
paper and straightforward if not explained otherwise. 



It proves convenient for the following discussions to 
transform the Hamiltonian into the quasi-particle basis, 
where it reads 



N 

EE 



(3) 



7!^' J denotes the annihilation [creation] operator of A G 
{-|-, — } quasi-particles in layer i e {a,p} with flavor / S 
{1,...,7V}. 

Since we are interested in conductivities we express the 
current in the particle-hole basis. As in the case of MLG 
the current decomposes into two contributions 



J = Ji Ji 



(4) 



II. MODEL & METHOD 

A. The model 

The model Hamiltonian in our problem in general con- 
sists of two copies of the free BLG Hamiltonian for the 

active and passive layer, respectively, as well as interac- 
tions within and in-between layers. It reads 



(1) 



a/p 



where a denotes the active layer and p the passive. 

denotes the free Hamiltonian in both layers, H^J^ the 
interaction within each layer, while H^^^ describes the in- 
teraction between layers. Disorder is implemented within 
each layer via iJ/nt- 

We restrict ourselves to a standard effective two-band 
description of BLG^^. Around the points K and K' we 
can expand the dispersion and end up with the following 
efiiective Hamiltonian: 



N 

Hi, = rf^k*'/W(k)*} ,with 



m = 0.054 me is the effective mass of BLG,^^ and mg is 

the electron mass. The two bands, subsequently identi- 
fled with electron (-|-) and hole (— ) band, have the respec- 

five energies, e^. — ji = — fi. In this approximation, 
the system is particle-hole symmetric at CN (/i = 0). / 
sums over the spin and valley degrees of freedom, con- 
sequently in BLG we have N = A. In the follwing, we 
maintain the parameter N as often as possible, to dis- 
tinguish it from other numerical prefactors. In the nu- 

iftl 

merical results it is always N = 4. ^fj" denotes the 
annihilation [creation] operator of layer i G {a,p} and 
flavor / e {!,••• , AT}. 



where the first accounts for the motion of the quasi- 
particles, while the second is the incoherent contribution. 
In terms of the electron and hole operators they read 



Ji 



e 

m 



E 

A/ 



E 



d^k 
(2^)2 



Ak7l/k)7;,/(k), 

(z X k) 



7;/k)7-/(k)-7i/(k)7+/(k) 



(5) 



where z is the unit vector perpendicular to the bilayer 
plane. Like in the case of MLG the d.c. conductivity is 
dominated by Jj corresponding to the quasi-particle con- 
tribution which produces a Drude peak in the absence of 
disorder and interactions. In the following we investigate 
the broadening of this quasi-particle peak. 



1. Coulomb interaction, screening, and disorder 

We consider the effect of Coulomb interaction, which 
in two dimensions is given by: 



V{q) 



eq 



(6) 



where e is the dielectric constant of the surrounding 
medium. The Coulomb interaction between the two 
layers separated by the distance d, which mediates the 
Coulomb drag, reads: 



eq 



-qd 



(7) 



Unlike MLG undoped BLG has a finite density of 

states at the CN point. Like in the two-dimensional elec- 
tron gas (2DEG), this is constant, and given by 



u{E) = N 



2m 



E \ =N- 



27r 



(8) 



3 



This implies that Coulomb interactions arc always 
screened and in the static limit Ii{q, to = 0) = v the 
Coulomb interaction within the random-phase approxi- 
mation (RPA) becomes 



VRPA{q) 



V 



1 + VU e{q + Nqtf) 



(9) 



where we have defined the inverse Thomas-Fermi screen- 
ing length 



(10) 



This corresponds to very small distances and defines the 
Bohr radius, which becomes important subsequently. 

The RPA for the two-layer setting is more involved 
and reads for the intra-layer interaction 



a/p 



v + {v^-u^)nj,,. 



2RPA 



(1 + viia) (1 + vup) - umjij, 



(11) 



' + {1- e-^i'^)NqTF 



e {q + NqTpf - e-'^i'^N'^q^p ' 



For the inter-layer interaction we find 

U 

U2RPA 



(1 + Vila) (1 + VILp) - um^iLp 



(12) 



qe 



e {q + NqppY - e-^9dN2q2^^ ■ 

In order to express the Coulomb interactions in the 
quasi-particle basis it proves useful to introduce the co- 
herence factor 

M,,v(k,q)=2(l + AA'^) , (13) 
where K = kx + iky and Q = qx + iqy Together with 

'7Ai,A2,A3,A4(k, p,q) 

= ^MA„A,(k,k + q)M;,„A3(p, P - q) (14) 



this allows to express the intra-layer interaction in the 
quasi-particle basis 



J2^ I (2^ I (^^•^i'^2'^3,A4(k,P,q)7l,A4,/(k + q)7l,A3,/'(P-q)7a,A2,/'(p)7a,Ai,/(k) . (15) 



For the inter-layer interaction we define 



^Ai,A2,A3,A4 (k, p,q) 



U{q) 



Mx,M (k, k + q)MA2,A3 (P, P " q) 



(16) 



where U depends on the distance d betweeen the layers. Then 
j2i, r j2„ r ^2 



rjap 



(2^ J (2^ J (2^^^i'^^'^^'^*^'''P'*l^ ^«.^4./^^ + 'l^^P.^3,/'(P-q)>.^2,/'(p)7a,Ai,/(k) . (17) 



If the system is at CN, there is no need to introduce im- 
purities in order to have well defined conductivities since 
the current carrying state is at zero total momentum and 
thus no momentum relaxation is needed. However, for fi- 
nite chemical potential the current carrying states excites 



a global momentum and we have to take into account 
breaking of translational invariance due to the presence 

of impurities to relax the momentum. We model disorder 
as charged impurity scattering with the Hamiltonian 



imp 



Ei:t.l^J 1^ .(|k - T+Nq.p) ^---^ q)7l,/('-)7v./(q) exp[^x. ■ (k, - k.)] . (18) 

i A, A' — i f — 1 



r 



III. KINETIC APPROACH 

We use the variational principle to calculate conduc- 
tivities from linearized Boltzmann equation as described 



in Chapter 7 of Ref. 26. Generalizations to drag are 
straightforward and can for instance be read in Ref. 24. 
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The one-particle distribution functions of electrons 
(+) and holes (— ) in the layer i G {a,p} is defined as 



/l(k,i) = (7i^/k,i)7l/(k,i) 



(19) 



We use the following ansatz for the non-equlibrium dis- 
tribution function 



1 



h=fx+ XX /a (1 - /a) -f^^ ■ eE , 



(20) 



where /° is the equilibrium distribution function (Fermi 
function) 



/A°(k) 



1 + e(^^-'*)/^ 



(21) 



We will use a one mode approximation, (20) , taking into 
account only the momentum mode ^ k. One can easily 
incorporate other modes like ~ k°,k'^,k^ into the for- 
malism. We have checked in the limiting cases of our 
study, that including these modes gives modified numer- 
ical values with corrections smaller than 1 per cent. Con- 
sequently, we have discarded them throughout. 

The collision integrals, which arc part of the Boltz- 
mann treatment, are shown explicitly in Appendix A. 
They can be derived directly from the interaction terms 
via Fermi's Golden rule or equivalently from a Keldysh 
calculation.^'' 



A. Drag resistivity 

We calculate the drag resistivity from the conductiv- 
ities obtained from a Boltzmann approach, sec App. A. 
We consider the response tensor which has a structure 
similar to the one in the Hall effect. The electric field E„ 
is only applied in the active layer o, and induces a current 
in the active layer called j^. It can also induce a current 
jp in the passive layer, which is the drag effect. Conse- 
quently, there are layer-diagonal and layer off-diagonal 
conductivities involved: 



<7d 0"p 







(22) 



This tensor includes the individual conductivities Ca and 
cTp of the active and of the passive layer, respectively. The 
layer off-diagonal part is called the drag conductivity, de- 
noted by (Td, and it corresponds to a cross-conductance. 
In the concrete experiment, however, the boundary con- 
ditions arc such that the passive layer docs not carry cur- 
rent. Demanding jp = requires a field counteracting the 



flow in the passive layer which is given by E„ 



iE 



This implies that the drag resistivity (or transresisitivity) 
is given by 



Like in the case of thermal transport pd can be fi- 
nite even if the individual conductivities Ua, <7p, and ad 
diverge. This is an effect of the boundary condition of 
vanishing charge flow in the passive layer analogous to a 
flnite thermal conductivity in thermal transport in Fermi 
Uquids.^^ 



IV. SINGLE LAYER CONDUCTIVITY 

In this section we describe the conductivity of an iso- 
lated BLG sheet. In the drag setup this limit is natu- 
rally achieved by setting the distance d between active 
and passive layer to inflnity. As in the case of MLG 
this quantity is well defined at the CN point, even in 
absence of impurities. This is a direct consequence of 
the particle-hole symmetry in the system, which implies 
that the current carrying state has effectively zero total 
momentum. 

The interaction limited single-layer conductivity has 
been investigated in MLG with the result^^ a = 

2 2 

0.76 -r-A^. In this expression a = corresponds 



?ia(T)2- — — -t- 

to the fine structure constant of graphene. It has been 
shown^^'^" that it has a logarithmic flow to zero upon 
decreasing temperature T. In the large N limit it was 
found^^ that the numerical prefactor is given by 0.69. 
While in MLG at CN it is justified to neglect the effect 
of screening due to the lack of density of states, this is 
not true for BLG. To see this one has to compare the 
temperature T against the energy associated with the 
Thomas-Fermi wavevector, which is given by 



Itf 
m 



xlO^K. 



(24) 



Thus for realistic temperatures one is always in the limit 
Tm <C q^p- Consequently, one always has to consider 
screening and we use the screened version of Coulomb in- 
teraction in the collision integral of the Boltzmann equa- 
tion. 

Solving the Boltzmann equation numerically in the one 
mode approximation yields the following expression for 
the minimal conductivity : 



27.4 34.0- 





27.4 0.259 — a/T [K] 



This implies that there is a constant conductivity in the 
limit of zero temperature and the leading temperature 
behavior is an unusual \/T. This can be rationalized as 
follows: assuming q^p ^ Tm the typical structure of 
integrals in the collision integral is given by. 



Pd = 



|EpP 

ja • Ep 



(23) 



dq 



qTpY 



5 
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FIG. 1: Interaction dominated single-layer conductivity. The 

leading temperature behavior is of the unusual \fT type. In 
the limit T — >■ 0, a converges to a fintic value a{T — >■ 0) = 

2 

27.4 The limiting value cannot be reached in experiment 
as a band gap opens below T w 5 K and our description 
breaks down.'^* 

where the exponential factor accounts for the Fermi dis- 
tribution function. The screened Coulomb interaction 
acts like a local interaction and dimensional analysis 
leads to a minimal conductivity which is independent of 
temperature, in agreement with the full numerical solu- 
tion. The leading correction is linear in q and conse- 
quently consequently produces a behavior corresponding 
to oc \fT (note that one can again understand this from 
dimensional analysis taking into account that the dynam- 
ical exponent is z = 2). There also is a high tempera- 
ture limit, in which the conductivity behaves linear in T. 
While this is a theoretical possibility this behavior dom- 
inates for temperatures above 10^ K, and consequently 
plays no role in experiments. Also, the zero tempera- 
ture limit cannot be reached in experiment since below 
T = 5 K BLG supposedly is gapped and our analysis 
breaks down. The result of the full numerical calculation 
of the conductivity is shown in Fig. 1, and the agreement 
with the \fT behavior is excellent within the relevant 
temperature range. 

V. COULOMB DRAG 

Motivated by recent experiments*^^""^^ many theo- 
retical studies of Coulomb drag in MLG have been 
published. ^^"^^'•^^"^^ We now investigate Coulomb drag 
in BLG and compare the results to those found for MLG. 

A. Qualitative results 

Our subsequent discussion has a number of natural di- 
mensionless parameters associated with different regimes. 

First and foremost, the parameter /i/T allows to distin- 
guish whether we are in the FL regime (/u/T ^ 1) or in 



the limit close to CN (/i/T <C 1). The parameter is 
always much larger than one for practical purposes, and 
consequently will not be discussed explicitly. We always 
assume Coulomb interaction to be screened. A further 
parameter important for the screening of the inter-layer 
interaction is given by dqTF- Due to qTF ^ 10^° m~^, 
which is the inverse Bohr radius, we are always in the 
limit dqxF ^ 1- The parameter which allows to distin- 
guish a clean system, meaning that the single-layer con- 
ductivities are dominated by Coulomb scattering, from a 
disordered system, where the single-layer conductivities 

are set by impurity scattering, is given by 5 = ^J^" . 
This corresponds to ratio of the impurity concentration, 
"imp, to the typical thermal momentum squared, or, 
equivalently, to the ratio of the scattering times asso- 
ciated with impurity scattering and Coulomb scattering. 

In order to get the overall qualitative picture we have 
plotted the full crossover curve of drag in Fig. 2. We 
observe that it is very reminiscent of the full crossover 
curve in MLG. However, as we will point out, there are 
important differences, mainly due to the fact that even 
in pristine BLG there is an additional energy scale due 
to the finite band mass. This mostly shows up in the 
distance dependence in the FL regime as well as on the 
behavior upon approaching the CN point. 




)U/T 

2 

FIG. 2: Drag resistivity for = 100 and dqrp = 1 (soUd 
line), dqTF = 2 (dashed line) , dqTF = 5 (dotted line). The 
single-layer conductivities Oa,<Jp are disorder dominated {g = 
1000). Then the maximum is located at /u ~ 2T. 



B. FL regime: /i/T > 1 

In the case of drag in MLG the effect of impurity scat- 
tering cancels exactly from pd in the FL regime. Con- 
sequently, all information about disorder is gone (this 
equally applies to a 2DEG). Using the same steps and 
manipulations as in the MLG,^^ one can show that within 
the one-band approximation 
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, irr. ,N ^4 1 1 

PdG«/T> 1 = 



duj 



(27r)2 sinh^ ^ 



I Uaj> (q, w)) |2 Im n++ (q, a;) Im n++ (-q, -u) , 



(26) 



where the density of states Ua^p was defined in Eq. (8) and 
Uap is the screened inter-layer interaction, Eq. (12). For 
1^1 ^ T the polarization function n++ (the +- 1- refers to 
the fact that only the quasi-particle contribution of the 
majority charge carrier is taken into account) assumes 
the standard Fermi-liquid form given by 



Imn++(g,w) 



kpq 



kpq 



- M 



m 



Scaling u coT and q ^ q/d allows to identify three 
natural dimensionless parameters: d = Tdm/kp, dkp 
and dqTF- We can then write the drag resistivity as a 
prefactor times a function which only depends on dimen- 
sionless parameters, namely 



h 4 



(27) with 



I kp 



F{d, dkp, dqTp) 



(28) 



2dkF 

F{d,dkp,dqTF) = J dq 



((g + NdqTp)^ - e-^iiNdqrpW J _ "" sinh" (w/2) 

-q/d 



q/d 

I 



dio- 



(29) 



r 



Since we have three independent parameters, we expect 
to find eight different regimes, depending on whether the 
parameters are much larger or smaller than unity. Two of 
these limits are internally incompatible with each other 
and there remain six regimes. Among these six only three 
are physical in the sense that for realistic systems they 
can be achieved. For more details regarding the analysis 
of the integral we again refer the reader to Ref. 24. 



1. d< 1 



while in the limit dqpp ^ 1 there is 
h 16n 



Pd = 



e2 3iV4 



{dqrpY 



(33) 



While the limit dqrp <C 1 theoretically makes sense it 
is not to be achieved in experiments due to the small 

screening length ^ q^]r- Consequently, there will be 
no saturation of distance dependence as in the case of 

MLG.24.39,40 



There are four limiting cases in this regime, two of 
them being physically sensible. We first discuss the sit- 
uation oi kpd ^ 1. In that situation we have either 
dqTP ^ 1 where the drag resistivity reads 

n 47r(, (6) II 1 ^ 

Pd 



1 



1 



{dqrpf {dkpf 



(30) 



or dqTP 1 where we find 

h 327r (tV (qrpV 



In 



1 



(31) 



\2NdqTP , 

If we go into the limit kpd <^ 1 we find in the limit 
dqrp < 1 

2 / \ 2 



h 327r fry (qTp\ 

/ NqpF + kp " 



X In 



TP 



Nq'_ 



TP 



kp 



(32) 



2. d> 1 

Here we find only two liming cases the reason being 
that this limit is incompatible with dkp <^1. As in the 
case of MLG, this regime is characterized by a linear in 
T behavior. In the limit kpd ;» 1 we have two situations, 
one in which dqTP ^ 1 with 



Pd = 



h 87r3 fT 



1 



1 



e2 15iV4 y y {dqTpY [dkpf 
while for dqTP <C 1 we have 



Pd = 



h32 fT\ fqTpV 1 



n \iJ, J \ kp J dkp 



(34) 



(35) 



Again it is important to note that the last regime is 
only theoretically possible. 
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3. Summary of results 

Prom the preceding discussion we conclude that in ex- 
periment only three different regimes can be realizable: 



Pd 



h_ 167r / T 
e2 3Afi I 



1 1 



h 47rC(3) t\ 

h 871-3 (t\ 1 ]_ 

e2 \n J {dqTFY (dkp)^ 



d <C 1 kpd <C 1 

d < 1 fc^d > 1 
d > 1 fci^d > 1 



(36) 



Unlike MLG there is no regime in which drag becomes 
independent of the distance d. 



C. Close to charge neutrality: /i/T <C 1 

In this limit drag resistivity is very sensitive to the 
ratio of the inelastic scattering time due to Coulomb in- 
teraction and the elastic impurity scattering time. This 
behavior has recently been analyzed for MLG^^'^^ and 
qualitatively the result in the BLG case is similar. The 
drag resistivity is well-defined and finite even in the clean 
system at finite chemical potential, where the individual 
layer and drag conductivities diverge. The parameter 
which allows to go from the dirty limit to the clean limit 
is given by 



frrii, 



'imp 



Tm 



(37) 



For g — > we realize the clean limit, while for g ^ oo 
we are in the dirty limit. In contrast to the FL regime, 
where elastic scattering drops out exactly, we find that 
the drag resistivity depends strongly on g. For disorder 
dominated single-layer conductivities, .g 3> 1, we find a 
maximum in the drag resistivity at ^ 2r, see Fig. 2. 
The same was shown to hold in MLG.^^'^^'^^ For smaller 
values of g, the maximum shifts to smaller values of n/T. 
In the following we use the parameter Jia,p = l^a,p/T 
for both layers, active and passive, individually as the 
dimensionless chemical potential. 

It was pointed out in the context of MLG^^'^^ that 
there are now two orders of limits in which one can de- 
scribe Coulomb drag in a clean system at particle-hole 
symmetry. If one first considers the system at CN, ^ = 0, 
drag vanishes irrespective of the disorder strength by 
virtue of lack of momentum transfer. However, if one 
first extrapolates to the clean system, — >• 0, at finite 
chemical potential and afterwards takes the limit // — >■ 
one ends up with a finite drag resistivity. Like in the 
case of MLG the first way to tak^ the; limits is the phys- 
ical one, resulting in a vanishing drag at particle-hole 
symmetry. In total, this shows that the drag resistivity 
sensitively depends upon the ration between inelastic and 
elastic scattering in the vicinity of CN. 

This behavior can best be exemplified in the unphysical 
zero distance limit, dqrp = 0, which we show in Fig. 3. 




0.1 0.2 0.3 0.4 0.5 0.67 



FIG. 3: Drag resistivity for a series of different values of g 

2 

in the idealized zero distance limit, dqxF = = 100). 

The single layer conductivities aa,crp range from interaction 
dominated g ^ 1 with a maximum at /i/T < 2 to disorder 
dominated for g ^ 1 in which case the maximum is located at 
H « 2T. For g — the ma:ximum pushes to zero. For /i > 2T 
the curves become independent of g and collapse to a single 
one, see section V B. 




0.1 0.2 



FIG. 4: Drag resistivity for = 100 and reasonable dis- 
tances dqTF = 10. The single layer conductivities aa,<^p here 
are mosty interaction dominated. This can be seen from the 

fact that the maximum is located at /i < 2T, in contast to the 
disorder dominated case, Fig. 2. For /x > 2T the curves be- 
come independent of g and collapse to a single one, see section 
VB. 



The general behavior assumes the following form: For 
small values of these parameters fXa and in the indi- 
vidual layers, we find a behavior which is well described 
by the function 



Pd « - 



' ai{g,dqTF, t^) + a2{g,dqTF, T^PaPi 



(38) 
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If we now keep g finite and consider \iia,p\/T <C 1 we find 



Pd 



h ai{g) 

MaMp hTm + b2nimp 



h T2 



(39) 



^imp 



where 61,62 depend on q^p/{Tm) and dqTF- This ex- 
pression goes to zero if we extrapolate the chemical po- 
tentials to zero. 

On the other hand, in the clean limit 5 — > we find 
ai 0. Now 



Pd 



P-aP-p 



h a2{g)PaPp 



e2 1 e2 



h 02(3) s->o h 



(40) 



meaning it extrapolates to a finite value as p — > 0. 
The constant c / 0, in contrast to the MLG^"*'^^, is not a 
universal number but it depends on q^p/{Tm) and dqTF, 
but not on g anymore (since we have performed the limit 
5 first). 

Wc have furthermore plotted a sequence of crossover 
curves with decreasing disorder strength in Fig. 4 for a 
more realistic value of dqrr = 10. As one can see form 
these curves, it is possible to extract information about 
both inelastic and elastic scattering times from Coulomb 
drag experiments carried out in the vicinity of CN. This 
has to be contrasted from the FL regime, in which infor- 
mation about disorder is absent. 



VI. CONCLUSION 

In this paper we have investigated interaction effects 
in transport phenomena in BLG. In a first part we stud- 
ied the minimal conductivity in pristine BLG. There wc 
find that the conductivity assumes a constant value in 
the limit T — >■ 0. Importantly, this result is valid in the 
clean limit and without lattice effects, meaning without 
Umklapp scattering. The first temperature correction 
to this is of the type VT. This has to be contrasted 



from the standard in Fermi liquids (neglecting ad- 

ditional logarithms and above all disorder) or the strictly 
independent of temperature (up to logarithms) behavior 
in MLG. In a second part we studied the Coulomb drag 
resistivity between two BLGs covering the whole range 
from deep within the Fermi liquid regime all the way to 
CN. We find that in the Fermi liquid regime drag behaves 
very similarly to drag in a standard 2DEG or MLG. In 
contrast to MLG, however, we find saturation of the drag 
resistance upon decreasing the distance d for physically 
reasonable parameters. In the vicinity of the CN point, 
like in the case of MLG, we find an interesting interplay 
between interaction effects and disorder making the drag 
resistivity strongly dependent upon the ratio of the two. 
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Appendix A: The collision integrals 

All conductivities are calculated using the variational 
principle. In Ref. 24 the formalism was adapted to cal- 
culate the drag resistivity. In the following wc present the 
collision integrals resulting from the Hamiltonians shown 
in section II. They are obtained by the linearization of 
the Boltzmann equation and taking into account all pos- 
sible second order processes.^'' 

The collision matrix has three contributions due to the 
three different scattering mechanisms: intralayer scatter- 
ing, interlayer scattering, and impurity scattering. We 
label it with the layer index a,p and particle-hole index 
A = ±, resulting in a 4 x 4 matrix. 

The layer diagonal part has contributions from all 
three processes 



aa,intra 



+ C 



aa^inter 
W 



+ SxyC 



(Al) 



The contribution from the elastic impurity scattering 
reads 



where M was defined in Eq. (13). For the intralayer Coulomb scattering we find 

Wa' mT2 7 (27r)2 7 (27r)2 7 (27r)2 ° i^2m 2m 2m 2m ^^^^^'^i'^^ 

X /^(fc)/°;"(Ai)(l - /A^"(|k + q|))(l - /";''(|ki - q|)) k. q ((1 - S^y) - <5.v) , (A3) 
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where 



R = 4Ar|T+__+(k, ki, q)|2 + |T+_+_(k, kj, ki - k - q)|2 - 4T+__+(k, ki, q)x 
X T|_+_ (k, ki, ki - k - q) - 4T|_+_ (k, ki, ki - k - q)T+__+(k, ki, q) , 



(A4) 



and T was defined in Eq. (14). Tlie contribution of the interlayer scattering to the layer diagonal part of the collsion 
matrix reads 



j^aa,inter 



2tt 

mr2 



J j2nrJ j2nrJ W)^ 



|k + q|^ 



2m 2m 2m 2m 
^n(k,ki,q) X /^(fc)/,r(fci)(l - /^(|k + q|))(l - /,V(|ki - q|)) k • q 

+ ^i2(k,ki,q)/°'"(fc)/,V(fci)(l - /°'P(|k + q|))(l - /°;«(|ki - q|))k. (k^AV - (ki - q)(l - 5xx')) 



+ 5 



2m 



kl 
2m 



|k + qp 
2m 



|ki 



2m 



-R2(k,ki,q) 



X /^(fc)/,V(fci)(l - /^(|k + q|))(l - /,V(|ki - q|))k.qMAv) 



where 



En=4Ar|T+__+(k,ki,q)|2, 
Ei2=4Ar|f+_+_(k,ki,ki-k-q)|2 , 
^2=4Ar|T++++(k,ki,q)|2, 

and T was defined in Eq. (16). Here only the large N processes contribute^*. 

The layer off-diagonal part has only a single contribution from interlayer scattering, which reads 



27r 



J 12^ J J2^J (2^ 



_ ^ _ |k + qp ^ |ki-qp 
\2m 2m 2m 2m 



i?n(k,ki,q) X /°'»(fc)/,V(fci)(l - /°'»(|k + q|))(l - flf{\k, - q|))k.q(l - Sxx') 



+ R,2{k, kl, q)/,^''^(fc)/,V(fci)(l - /°^P(|k + q|))(l - f?f{\k, - q|)) k. (ki(l - d^x') - (k + q)5,v) 



+ s(— + 



,2 II, I „\2 



I k'i |k + q| 
V2m 2m 



|ki 



i?2(k, ki,q) 



2m 2m 

X /,^''^(fc)/,r (A:i)(l - f^'^ilk + q|))(l - /,r (|ki - q|)) k • q (S.y) 



(A5) 



(A6) 
(A7) 
(A8) 



(A9) 



To calculate the conductivities using the variational prin- For clarity we have omitted the particle-hole index, 
ciple, we need the driving terms of the relevant modes^^. 

For the model at hand and the momentum mode, that 
we have used, they are of the form 

= (01,01) , and DP = {Dl,DP_) , (AlO) 

where 

^±'=1^ / (0/±"''(l-/±"'')'^-^ • (All) 
Then the conductivities can be calculated according to 



(jaa (jap \ "1 / jy 



'^'^ ~ h T ' ycp'' cpp ) ' \ Q 

^'^ rr. r^v^ ( C^P f 

= (0, UP) ■ 
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